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This paper studies the mutation-selection balance in three simplified replication models. The first 
model considers a population of organisms replicating via the production of asexual spores. The 
second model considers a sexually replicating population that produces identical gametes. The third 
model considers a sexually replicating population that produces distinct sperm and egg gametes. All 
models assume diploid organisms whose genomes consist of two chromosomes, each of which is taken 
to be functional if equal to some master sequence, and defective otherwise. In the asexual population, 
the asexual diploid spores develop directly into adult organisms. In the sexual populations, the 
haploid gametes enter a haploid pool, where they may fuse with other haploids. The resulting 
immature diploid organisms then proceed to develop into mature organisms. Based on an analysis 
of all three models, we find that, as organism size increases, a sexually replicating population can 
only outcompete an asexually replicating population if the adult organisms produce distinct sperm 
and egg gametes. A sexual replication strategy that is based on the production of large numbers 
of sperm cells to fertilize a small number of eggs is found to be necessary in order to maintain 
a sufficiently low cost for sex for the strategy to be selected for over a purely asexual strategy. 
We discuss the usefulness of this model in understanding the evolution and maintenance of sexual 
replication as the preferred replication strategy in complex, multicellular organisms. 
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I. INTRODUCTION 



A. Why sex? 



The preference for sex as the exclusive mode of replica- 
tion in complex multicellular life is a long-standing prob- 
lem in evolutionary biology P, 0, H, 0, The difficulty 
with understanding the preference for sexual replication 
is that, at first glance, sex appears to be a cumbersome 
and inefficient mode of reproduction: While asexually 
replicating organisms can simply produce asexual spores 
that then develop into new adult organisms, sexual repli- 
cation involves the mixing of genetic material from two 
distinct organisms. 

The need to combine genetic material from two dis- 
tinct organisms incurs a fitness penalty, in the form of 
time and energy costs that each sexually replicating or- 
ganism must pay in order to find a genetic recombination 
partner. Further, in cases where a population employs 
a sexual replication strategy that involves distinct sexes 
producing distinct egg and sperm gametes, the potential 
rate of reproduction is half that of an asexually replicat- 
ing population, where each organism produces diploid 
eggs. This 50% difference in reproduction rates is known 
as the two-fold cost for sex . 
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B. Theories for the existence of sex 

Despite the disadvantages for sex, its preference, which 
is in many cases exclusive, in the higher organismal lines 
indicates that it must confer an overall fitness advantage 
in complex organisms. 

A variety of theories have been proposed for the bene- 
fits of sex, which generally fall into one of two categories: 
(1) Sex increases adaptability; and (2) Sex prevents the 
accumulation of deleterious mutations. 

The first category of theories, originally due to Weis- 
mann, state that sex evolved because it allows small pop- 
ulations to adapt more quickly to changing environments 
0. By allowing for recombination amongst different or- 
ganisms, sex potentially allows for isolated beneficial mu- 
tations to become incorporated into a single organism. 
The result is that sex can potentially greatly speed up 
rates of adaptation. For complex, slowly replicating or- 
ganisms, or in small populations, this feature of sex can 
provide a significant fitness advantage. 

The second category of theories, by contrast, holds 
that sex evolved because it prevents the accumulation 
of deleterious mutations. According to this theory, 
sex allows organisms to discard defective genes in their 
own genomes and replace them with functional copies 

0,1 IE, [Hi [a. 

Within each category there are a number of compet- 
ing theories that are variations of a general theme. In 
the context of the first category, the two most common 
theories are the Vicar of Bray Hypothesis^ and the Red 
Queen Hypothesis. 

The Vicar of Bray Hypothesis simply assumes that sex 
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allows for faster adaptation in dynamic environments. 
The name is derived from an English cleric who would 
change his political and religious views as necessary in 
order to remain in office [l| . 

The Red Queen Hypothesis is a somewhat more com- 
plex version of the Vicar of Bray hypothesis. It states 
that sex provides a fitness advantage as a result of a 
constant co-evolutionary genetic "arms race" with fast 
replicating and evolving parasitic organisms [1, Q . The 
name is derived from a character, the Red Queen, in the 
story Through the Looking Glass by Lewis Carroll. In the 
story, the Red Queen states, "It takes all the running you 
can do, to keep in the same place." Q 

In the context of the second category, the two most 
common theories are the Genetic Repair Theory, and the 
MuUer's Ratchet Theory. The Genetic Repair Theory 
simply states that sex prevents the accumulation of dele- 
terious mutations 0, B [ll- The MuUer's Ratchet The- 
ory argues that sex slows down a mutation-accumulation 
phenomenon in small populations known as MuUer's 
Ratchet (see Figure 1). 

It should be noted that the various theories for sex need 
not be contradictory. That is, sex may indeed prevent 
the accumulation of deleterious mutations and allow for 
faster adaptation in dynamic environments. The reason 
for this is that the prevention of accumulation of delete- 
rious mutations is accomplished by discarding defective 
genes, while adaptation occurs by bringing together ben- 
eficial genes. Since these two processes are similar, if not 
exactly equivalent, it is possible that both categories of 
theories may provide complementary, rather than rival, 
explanations for the existence of sex. 

Furthermore, although a small population may not be 
a requirement for a selective advantage for sex, the se- 
lective advantage for sex may nevertheless be greater in 
small populations than in larger ones. 



C. Problems with the various theories for sex 

The four theories for the selective advantage for sex 
have certain difficulties that makes each of them incom- 
plete. The adaptability category of theories requires a 
dynamic environment for sex to confer a fitness advan- 
tage. However, a number of sexually replicating organ- 
isms (sharks and crocodiles, for instance), have remained 
apparently unevolved for tens of millions of years in 
what are seemingly fairly static environments jl3l.[T3.[lH. 
While there may be some environmental dynamics that 
is difficult to detect (parasites, for instance), it is not 
immediately clear that a dynamic environment is a nec- 
essary condition for the emergence and maintenance of 
sex. 

The MuUer's Ratchet Theory suffers from its reliance 
on a small population. This is an ill-defined term, since in 
one context a given population may be considered large, 
while in another context it may be considered small. For 
example, is the current human population of approxi- 



mately 7 X 10 people considered small or large? If this 
is a large population, then the MuUer's Ratchet Theory 
would argue that the human population should eventu- 
ally become asexual. While this supposedly large human 
population has existed for much too short a time com- 
pared with macroevolutionary time scales, it is neverthe- 
less unclear if the human population would eventually 
become asexual if maintained at current population size 
and density. 

The Genetic Repair Theory is the generally accepted 
theory for the existence of sex, mainly because it is the 
one that requires the fewest assumptions. Nevertheless, it 
too is problematic, for if sex prevents the accumulation 
of deleterious mutations, then why don't all organisms 
replicate exclusively sexually [T6|? 

The issue here is not why genetic recombination is ad- 
vantageous. It is well-known that genetic recombination 
between organisms occurs at all levels of organismal com- 
plexity, including bacteria and even viruses [ITj]. Fur- 
thermore, certain single-celled organisms, such as Sac- 
charomyces cerevisiae, or Baker's yeast, can engage in 
a form of sexual replication as part of a stress response 
[isj . However, the single-celled organisms that are capa- 
ble of genetic recombination either replicate themselves 
asexually, or only replicate sexually when stressed. This 
behavior is in sharp contrast to more complex, multicellu- 
lar organisms, which replicate sexually either exclusively 
or nearly exclusively. 

Clearly then, there are regimes where the advantages 
for sex are outweighed by its disadvantages. A complete 
theory for the existence of sex must be able to identify 
the regimes where either sexual or asexual replication 
are respectively dominant, in a manner that is consistent 
with observation. 



D. Haploid fusion time and its affect on the 
selective advantage for sex 

1. A density-independent model assuming first-order 
haploid fusion kinetics 

In [l^ , the author developed a mathematical model de- 
scribing the evolutionary dynamics of a unicellular pop- 
ulation that could replicate either asexually or sexually. 
The sexual replication mechanism was loosely based on 
the sexual stress response in Saccharomyces cerevisiae 
(Baker's yeast) (see Figure 2). The sexual replication 
model explicitly incorporated a time cost for sex, in the 
form of a characteristic time for haploid cells to find each 
other and recombine. 

What the author found was that when the haploid 
fusion time was small compared to the characteristic 
growth and doubling time of the cell, sexual replication 
was the advantageous strategy. However, when the hap- 
loid fusion time was large compared to the characteristic 
growth and doubling time of the cell, then sexual repli- 
cation was the disadvantageous strategy. These results 
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implied that sex should be the preferred replication strat- 
egy in slowly replicating organisms, and, all other fac- 
tors being equal, in regimes of high population density. 
These results therefore delineated regimes where sexual 
and asexual replication are respectively advantageous in 
a biologically consistent manner. 



2. A density- dependent model assuming second-order 
haploid fusion kinetics 

The sexual replication model developed in [l^ as- 
sumed a density-independent haploid fusion time, an 
unrealistic simplifying assumption that was changed in 
[20| to a second-order rate process (the assumption of 
a second-order rate process is a natural one to make if 
haploid fusion is modeled as occurring via random bi- 
nary collisions between haploid cells). While the results 
of this paper were qualitatively similar to the results 
in [l^, the paper raised interesting questions regarding 
the emergence of sexual replication in an asexually repli- 
cating population that did not arise with the density- 
independent model. 



3. Multicellular organisms 

Although the work in [l^[23| suggested that sex should 
be favored in slowly replicating organisms, it did not work 
with a replication model appropriate for multicellular or- 
ganisms. With multicellular organisms replicating asex- 
ually, single-celled spores, or eggs, are continually pro- 
duced by the adult organism, which may then develop 
into an adult organism. 

In the case of multicellular organisms replicating sexu- 
ally, haploid gametes are continually produced by diploid 
germ cells. These gametes may recombine with other 
haploid gametes, to form "post-fusion" spores (e.g. fertil- 
ized eggs) that may then develop into an adult organism 
(see Figure 3). 

In [2l|, the author developed mathematical models 
for asexual and sexual replication based on sporulation. 
The sexual model did not assume internal fertilization 
or even explicit mating between organisms. Rather, the 
model adopted a sexual replication model based on ex- 
ternal fertilization, whereby haploids are released into 
an external aqueous environment, and fuse with other 
haploids in the haploid pool. More specifically, the sex- 
ual replication model assumed a specific form of external 
fertilization known as broadcast fertilization or spawning. 
In this fertilization mechanism, the adult organisms do 
not attempt to get closer to one another before releasing 
their gametes. Rather, the adult organisms simply re- 
lease their gametes into a general pool, so that the adult 
organisms do not facilitate haploid fusion in any way [2^ . 
This model reflects the fact that terrestrial life, both uni- 
and multi-cellular, developed in the oceans, and that sex- 
ual replication in the older organismal lines involves the 



release of haploids into the water Therefore, the 

simplest models capable of identifying a selective advan- 
tage for sexual replication in complex organisms, and a 
disadvantage for sexual replication in simple organisms, 
should be able to do so with this primitive fertilization 
mechanism. 

Assuming second-order kinetics for the haploid fusion 
rate, the author found that sex is indeed the preferred 
replication strategy when the time cost for sex is small 
compared to the characteristic growth time of the organ- 
ism. However, the model also suggested that as organ- 
ism size increases, the consequent decrease in population 
density implies that the cost for sex should increase as 
organism size grows. It was found that, in this case, sex 
can only provide a selective advantage over asexual repli- 
cation if the rate of production of gametes decreases as 
1/iV, where N is the number of cells in the adult organ- 
ism. This is clearly a problematic result, since multicellu- 
lar organisms can produce gametes in great numbers (the 
males of some species, for instance, can produce millions 
of sperm every day) [23[. 



E. Sexual replication with sperm and egg gametes 

As will be shown in this paper, it turns out that the 
resolution to this discrepancy rests on the following cru- 
cial assumption: If a fertilized spore (egg) develops in an 
external aqueous environment, then it must contain a suf- 
ficient amount of material in order to allow the initial cell 
to develop into an immature multicellular organism that 
is capable of functioning independently of its parents. In 
short, the adult organisms do not produce gametes that 
are of microscopic size, but rather eggs that contain a 
sufficient amount of extra-genetic material (the "yolk" ) 
to allow the organism to develop sufficiently so that it 
has a reasonable chance of survival after hatching. 

Working from this assumption, we find that, as organ- 
ism size increases, it is necessary for the adult popula- 
tion to produce two distinct types of gametes: Large, 
relatively immobile eggs that contain the necessary non- 
genetic material to produce a viable immature organism, 
and small, highly mobile sperm that fertilize the eggs. 
This type of replication strategy is in contrast to the 
replication strategies considered in [l^, [2l[ , which as- 
sumed gametes of more or less identical size (isogamy) . 

Distinct sperm and egg gametes (anisogamy) are found 
to be crucial for maintaining a sufficiently low cost for sex 
so that it is the preferred replication strategy as organism 
size increases. The results of this distinct gamete model 
therefore provide a useful insight into the preference for a 
sexual replication strategy in complex multicellular life, 
and also suggests a resolution to the two-fold cost for sex 
and the necessity of males. 

This paper is organized as follows: In the following 
section (Section II), we review the asexual replication 
model considered in plj . Since we are explicitly con- 
sidering eggs in this paper, some of the terminology is 
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slightly changed. In Section III, we review the identical- 
gamete sexual replication model considered in 21|. In 
Section IV, we develop and analyze the distinct-gamete 
sexual replication model. In Section V we compare all 
three replication models, and show that a distinct-gamete 
model is necessary for sexual replication to be advanta- 
geous over asexual replication as organism size increases. 
In Section VI we further discuss the implications of our 
results. Finally, in Section VII we summarize the main 
conclusions of this paper, and describe plans for future 
research. 

As a concluding note for this section, it should be noted 
that there has been a considerable amount of research on 
the selective advantage for a distinct-gamete, or anisoga- 
mous, sexual replication strategy 12, l^,!!!,!!^. Indeed, 
one of the main theories for the selective advantage of an 
anisogamous strategy is that it reduces the time cost as- 
sociated with haploid fusion. 

In one version of this theory, an anisogamous strategy 
reduces the time cost associated with haploid fusion be- 
cause many sperm can be produced for every egg, so that 
the eggs are essentially bathed in a "sperm cloud" [23 |. 

In another version of this theory, an anisogamous strat- 
egy reduces the time cost associated with haploid fusion 
because the small sperm are highly mobile and the eggs 
provide a large surface area for contact, which together 
leads to an increase in sperm-egg contact frequency [26l |. 

However, previous research on the evolution of 
anisogamy has focused on the selective advantage of a 
distinct-gamete sexual replication strategy over identical- 
gamete, or isogamous, sexual replication strategies. This 
paper, by contrast, argues that, within the context of 
broadcast fertilization, gamete differentiation is neces- 
sary in order to provide a selective advantage for the 
sexual strategy itself. Nevertheless, because the idea 
that distinct gametes reduce the time cost for sex has 
been advanced before, our arguments will be similar to 
those of previous authors. In particular, this paper fol- 
lows the approach of Dusenbery, and applies the kinetic 
theory of gase s to analyze the haploid fusion process 
[2i,[23,[2i,l2i,[3|. 



II. THE ASEXUAL REPLICATION MODEL 



A. Definitions 

We begin by summarizing the asexual replication 
model considered in 12 1| . and which will be used as a 



basis for comparison with the other replication models in 
this paper. 

We assume that we have a population of diploid organ- 
isms, whose genomes consist of two chromosomes. Each 
chromosome may be represented as a linear symbol se- 
quence a = siS2 ■ ■ ■ sl, where L is the number of letters 
(equivalently, bases) in the sequence, and where Si de- 
notes the i*^ letter (so Si = A, T, G, or C for DNA, A, 



U, G, C for RNA). A given genome may then be repre- 
sented by the set {a, tr'}, where a and a' denote the base 
sequences of the two chromosomes. 

We also assume that there is a "master" sequence, de- 
noted (To 5 for which the chromosome is functional, and 
that any mutation to ctq renders the chromosome defec- 
tive. It is then assumed that the fitness of an organism 
(to be defined below) is determined by the number of 
functional chromosomes (zero, one, or two) in the organ- 
ism. 

While this simplifying assumption is clearly a highly 
coarse approximation of organismal fitness, it is the gen- 
eralization of the Single Fitness Peak Landscape com- 
monly used in quasispecies theory [3l|, [13, [H, |3^ . [35| . 
This fitness landscape is used because it is the simplest 
landscape that is both analytically tractable, and that 
gives phenomenological results that are in qualitative 
agreement with a number of important biological effects 
(in certain cases, these results can even be quantitative) 
[3a |. Therefore, we regard this choice of landscape as a 
useful starting point for our model. It should be men- 
tioned, however, that recent research in quasispecies the- 
ory [33, [H, [3^, and evolutionary dynamics in general, 
has focused on developing and analyzing more realistic 
fitness functions [13, EH . 

Within the context of our chosen fitness landscape, 
the population can be divided into three categories: The 
first category consists of organisms with genome {ao, <Jo}- 
These organisms have genomes consisting of two func- 
tional chromosomes. If we define the chromosome ao to 
be viable, then we may define organisms with two viable 
(i.e. functional) chromosomes to be of type vv. 

The second category consists of organsims with genome 
{(To, cr ^ Co}- These organisms have genomes consisting 
of one functional and one non-functional chromosome. 
If we define the chromosomes cr ^ (Tq to be unviable, 
then we may define organisms with one viable and one 
unviable chromosomes to be of type vu. 

The third and final category consists of organisms 
with genome {a ^ iTo,(t' 7^ ctq}. These organisms have 
genomes consisting of two non-functional chromosomes, 
and are defined to be of type uu. 

The asexual replication process, as illustrated in Fig- 
ure 4, is assumed to occur as follows: An adult organ- 
ism continually generates asexual, single-celled spores 
at some genome-dependent rate uj^g.,jiy. These spores, 
which we also refer to as immature organisms, develop 
into adult organisms with a first-order rate constant, de- 
noted K{a,a'}- 

We define ujyy = w^^^^^^}, uj^u = i^{ao^a^,7o}> and 
w„„ = uj{a^c,o.a'^cro}- Wc also define k„„, and Kuu 
in an analogous manner. 

Throughout this paper, we will assume, unless other- 
wise stated, that, = lOuu < ^^vu < ^w, and — k„„ < 
Kvu ^ Kvv This makes sense, since organisms with two 
defective chromosomes should not be expected to grow 
or sporulate, and organisms with one defective chromo- 
some should not be expected to grow or sporulate more 
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rapidly than organisms with no defective chromosomes. 



Furthermore, we define a = Kvu/kw, and (3 = 
ujyu/'-^vv These parameters may be interpreted as the 
growth and sporulation penalties, respectively, associ- 
ated with having a defective chromosome. 

The actual production of spores is assumed to occur as 
follows: A subset of the cells in the adult organism are 
responsible for producing spores. These are the spore 
stem-cells, otherwise known as germ cells. Each of these 
cells periodically divide to produce two daughter cells, 
one of which remains with the population of the spore- 
producing cells, while the other goes on to become a 
spore. 

Within this model for spore production, we assume 
that a given parent chromosome has a probability p of 
replicating correctly. We also assume that the sequence 
length L is sufficiently large that the probability that an 
already mutated base will mutate again is negligible, so 
that any new mutations in the daughter strand must oc- 
cur in a previously un-mutated portion of the genome. 
This assumption, termed the neglect of backmutations, 
means that the probability of a w parent chromosome 
producing a v daughter chromosome is p, the probabil- 
ity of a w parent chromosome producing a u daughter 
chromosome is 1 — p, and the probability of a u parent 
chromosome producing a u daughter chromsome is 1. 

Finally, we assume that when a given spore-producing 
cell divides, the daughter cell that remains the spore- 
producing cell retains the parent chromosomes, while 
the daughter cell destined to become a spore retains 
the daughter chromosomes. This chromosome segrega- 
tion mechanism, termed immortal strand segregation, is 
known to occur in adult stem cells, and is believed to 
be the chromosome segregation method in budding yeast 
[12, m, S) mi • Because of its apparent ubiquity in stem 
cell populations, we adopt this chromosome segregation 
method here, though it is possible to consider random 
recombination, as was done in [l^ [20l|. 

Figure 5 illustrates the details of the spore-production 
process. 



B. Population genetics equations and the mean 
fitness 



If we let Ua 



denote the number 



of adult organisms with genomes vv, vu, and uu respec- 
tively, and if we let riai^w, nai^vu, nai^uu denote the num- 
ber of immature organisms (spores) with genomes vv, vu, 
uu respectively, then the time evolution of the various 
populations is given by the following system of ordinary 
differential equations: 



drin 



dt 



dn„ 



dt 



dn„ 



dt 



^VV I ''ai.VV 



^vu I ''ai.vu 



dn 



ai,vv _ 2 



dUn 



dt 



dUn 



dt 



2j>(l P)^vv'^am,vv P^vu'^am,vu ^vu^ai, 



(1) 



The central object of interest, that we will obtain 
from the equations, is the mean fitness of the popu- 
lation, which is time-varying and therefore denoted by 
ii{t). If we let n = Uam.w + nam.vu + nam,uu denote 
the total population of adult organisms, then we define 
ii{t) = (1 / n) {dn / dt) , so that the mean fitness is simply 
the per-capita growth rate of the population as a whole. 

The mean fitness measures the effective first-order 
growth rate constant of the population. As a result, if two 
separate populations are placed in a given environment, 
the population with the higher mean fitness will drive the 
other to extinction. Therefore, the mean fitness allows 
us to determine which replication strategy (e.g. asex- 
ual, sexual with identical gametes, sexual with distinct 
gametes) is dominant in a given parameter regime. 

It should be noted that we are working with a group se- 
lection approach for determining which replication strat- 
egy is dominant for a given set of parameters. While 
such an approach is often the easiest to analyze, one has 
to be careful in using it, since it is well-known from evolu- 
tionary game theory that the evolutionarily stable strat- 
egy (ESS) does not always coincide with the replication 
strategy that maximizes the overall fitness of the entire 
population [2^ [4^ . This is an issue that we will discuss 
toward the end of this paper. 

Although in this simplified model the population can 
grow indefinitely, the population fractions eventually 
reach a steady-state, and so the mean fitness converges 
to a steady-state as well. To determine the steady-state 
mean fitness, we define a new set of parameters that con- 
verge to steady-state as well, and re-express our dynam- 
ical equations in terms of them. 

To this end, we define Xaq.rs = naq^rs/n, where q = 
m,i, and rs — vv,vu,uu. With these definitions, we 
have that. 



R{t) 



^uw^am.uu 



(2) 



Re-expressing the dynamical equations in terms of the 
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Xaq.rs population ratios, we obtain, 

dXam.vv 



dt 

dXam.vu 



dt 



dx„ 



dt 



dx. 



P ^vvXam,vv \f^vv \ f^\J^))Xai.vv 
2j>(l — p)^vvXam,vv ~l~ P^vuXarn^vu 



dt 



dt 



dxgj^uu _ _ \2 M (^ — ^ 

J, — P} ^vvXam,vv i l-L P)^vuXam,vu 

at 



C. steady-state behavior 

The steady-state solution may be found by setting the 
left-hand-sides of these equations to zero, and solving 
for the various population ratios and for the steady-state 
mean fitness. We obtain [2ll |. 

R{t — oo) 

•^ai.vv — Xam.vv 

R{t = oo) 

Xai,vu Xam,vu 



K,{t — oo) 



(4) 



which may be substituted into the bottom three equa- 
tions to give, 

= \l'^{t = oo) + Kyt^K{t = oo) p (jJvvl^vv]Xam,vv 
= -2p(l - p)uJvvKvuXam,vv + [^(t = Oo)^ 

^vu^i^t oo) P^vu^vu\Xam,vu 
— (1 p) ^vv^uuXajn.vv (1 p)^vu^uuXam.vu 
-\-[K(t — Oo)^ -|- Kuu^{t = oo) — <JJuufiuu\Xam,uu (5) 

We have three cases, giving rise to three distinct 
steady-states: 

1. Case 1: Xam,vv > 



If Xam,vv = at steady-state, then the second-to-last 
equation reads, 

= {fi{t = OO)^ + KyuK{t = oo) - pUJvuKvu)Xam,vu (8) 

If Xam.vv = and Xam.vu > at steady-statc, then this 
implies that. 



K{t = oo) = Ka.2 = + a/^ + 4— p) (9) 

2 y Kv 



3. Case 3. X am.vv — •^am.vu — 

If Xam.vv and Xam.vu = at stcady-state, then the last 
equation reads, 

= {R{t = OO)^ + K„„K(i = oo) - UJuuKuu)Xam,uu (10) 



Since Xa 



1, then x„ 



^am.vv ~1~ Xam.vu ~\~ Xam.uu 

Xam,,vu = ^ Xam.uu = 1, and hcncc the steady-state 
mean fitness must be given by. 



Hit - oo) - Ka.3 = ^Ku«(-1 + \/l+4 — ) (11) 



Based on a stability analysis of the steady-states, it is 
possible to show that R(t = oo) = max{Ka.i, Ka.2, f^a.a}- 
Using the standing assumption that ujuu — = 0, then 
Ka.3 — and R{t ~ oo) = max{Ko_i, Ka.2}- 

In Appendix A, we show that there exists a Pcrit £ 
[a(3,(3] such that the mean fitness is given by Kq^i for 
p > Pcrit, and by Ka^2 for p < Pcnt [SI. 

The transition at Pcrit is associated with a localiza- 
tion to delocalization transition over one of the chro- 
mosomes in each organismal genome, so that mutation- 
accumulation in one of the chromosomes is governed by 
random genetic drift. This phenomenon is well-known in 
quasispecies theory, and is termed the error catastrophe 



III. SEXUAL REPLICATION WITH 
IDENTICAL GAMETES 



If Xam.vv > at steady-state, then it is possible to see 
that the steady-state mean fitness must satisfy, 

— K{t — OO)^ + KyyR{t — oo) — p^OJvvKvv (6) 

so that, 

R{t = oo) = Ka.l = \-Kyy{-l + 1/1+ 4— p2) (7) 

Z y K'ny 



We now turn our attention to the identical-gamete sex- 
ual replication model analyzed in J21j. Here, we assume 
that when an initial diploid spore is formed from the di- 
vision of a spore stem-cell, it splits into two haploid ga- 
metes that then enter a haploid pool and fuse with other 
haploids. The resulting diploid is a "post-fusion" spore 
(i.e. a fertilized egg) that we also term a new imma- 
ture organism, which then develops into a mature adult 
organism in a similar manner to the asexual organisms. 
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A. The basic model 

As with the asexual rephcation model, we have sporu- 
lation and maturation rate constants denoted by w^^, 
ujyu, ^uu, and Kyy, Kyui i^uw Because there is no reason 
to assume otherwise, we will assume that these values are 
identical to the corresponding values in the asexual pop- 
ulation (the populations are taken to be identical in all 
respects except in one aspect of their replication cycles). 

We model haploid fusion as a bimolecular collision re- 
action characterized by second-order kinetics, with a rate 
constant denoted 7. We also assume that only the v 
haploids can fuse with one another, since the u haploids 
have a single defective chromosome, and therefore they 
can presumably no longer take part in the replication 
process. This restrictive assumption will be relaxed in 
future research when we consider other kinds of mating 
strategies (such as random mating). 

Since only the v haploids can fuse with one another, 
this replication strategy implies that only vv genomes are 
present in the population. We therefore need only con- 
sider the population of mature vv organisms, immature 
vv organisms, and v haploids. If we let the numbers of 
these respective populations be denoted by nsm.w, nsi^vv, 
and riu, respectively (the "s" stands for sexual), then, if 
V denotes the system volume, the evolutionary dynamics 
of the population is governed by the following system of 
ordinary differential equations: 



expressing the evolutionary dynamics equations in terms 
of these parameters, we obtain. 



dt 



^ = 2W„„P - JPXI - K{t)Xy 



(13) 



It should be noted that this version of the identical- 
gamete model ignores gamete death. 



B. Steady-state mean fitness 

Following [2l[, we may find the steady-state mean fit- 
ness for this model by setting the left-hand sides of the 
above two equations to 0. We obtain. 



12'^'"' {1 I ^'■^ ^ oo) ^j ^ft = 00) ^^^^ 



IP 



which may be substituted into the second equation to 
give, after some manipulation, 



(iiVsip — (1 -I- '^C*^"") ^ K{t=oc) -^2 



1 



_ (15) 
^ K{t=co) _|_ K(t=oo)-j 2 7p 



dt 



dn 



dt 



17^ 



si,vv ^ f ^ \ 2 



dn. 



dt V 



(12) 



where the factor of 1/2 in the second equation comes from 
the fact that two haploids form one immature organism, 
so that if (7/1^)71^ is the rate at which haploids disappear 
due to fusion, then (l/2)(7/T^)ri^ is the rate at which 
immature organisms are produced. 

Furthermore, the factor of 2 in the third equations 
comes from the fact that the diploid spores divide in 
two, so that the rate of production of haploid gametes 
is twice the rate of production of diploid spores. Since 
a given chromosome is replicated correctly with a proba- 
bility p, the rate of production of viable haploids is given 
by p times the rate of production of haploids. 

Now, we make the further assumption that the adult 
organisms take up a certain amount of space, so that, as 
the adult population grows, so does the system volume. 
To this end, we define a population density of adult or- 
ganisms p ^ ns„i,vv/V, and assume that p is fixed. Fur- 
thermore, we define the population ratios Xsi.w and x^ 

via Xsi^vv — '^sijVV /^sm^vv : and Xy — / Tlgjji^yy . 

Since the adult population consists only of vv organ- 
isms, it follows that n = Usn 



ness K{t) = {l/usm /dt) 



and so the mean fit- 
Re- 



A closed form expression for the mean fitness may be 
found in the limit Kyy/{'yyyp) — s- 0. This corresponds 
to a situation where there is no cost for sex, since the 
characteristic growth time, 1/Kyy, is very large compared 
to the characteristic haploid fusion time, 1 / (7p) . 

In this case, the mean fitness is obtained by solving the 
quadratic. 



^ R{t^ 00) ^^ ^ R{t = 00) UJyy^_^ ^^g^ 



SO that. 



k(< = oo) = -k™(-1 + J1 + 4p— ) (17) 



KyyX 



which may be shown to be larger than 1 for p G (0, 1), 
and larger than Ka.2 as long as p > and either a < 1 or 
/3 < 1. 

Therefore, when there is no cost for sex, it is the pre- 
ferred replication strategy. 



C. Asexual versus sexual replication as a function 

of K,yy/{jp) 

We may show that K{t — 00) /Kyy is a decreasing 
function of Kyy/"fp. To see this, note that the function 
f{x,X) = {X— {l + x)x)'^ /{x^{l + x)) is a decreasing func- 
tion of a; as a; increases from to the value of x, denoted 
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X* , for which = A — (1 + x)x. As x increases from to 
X* , the function f{x,X) decreases from oo to 0. 

Therefore, as K^v/iip) increases from to oo, the mean 
fitness K{t = oo) decreases from its value when there is no 
cost for sex down to 0. Because the mean fitness of the 
sexual population is larger than the mean fitness of the 
asexual population when there is no cost for sex, it fol- 
lows that there is a unique value for K««/(7p), denoted 
{K'vv/{ip))crit, where the asexual and sexual mean fit- 
nesses become equal. Below this value of k,vv/{ip)i the 
sexual mean fitness is larger, while above this value of 
Kw/iip), the asexual mean fitness is larger. 

To determine the critical value of n^v/iip) where the 
asexual and sexual strategies yield identical mean fit- 
nesses, we simply substitute the asexual mean fitness into 
Eq. (15), and this directly gives us our result. 

We can study the asymptotic behavior of 
(k^v / il p)) crit in both the small and large uj^y/K^y 
limits. When ujyy/Kyy is small, then the asexual mean 
fitness is given by, 



K{t = oo) IjJyy 



-XP 



(18) 



where x = P f^^ depending on whether the asexual 
mean fitness is given by Ka,i or Ka^2- 

In this case, keeping the lowest-order terms involving 
i^vv/nyy in the numerator and denominator of Eq. (15), 
we have. 



,Kyy ^) 

I crit — ^ ^ 

IP X P l^vv 



(19) 



When ujy^l Kyy is large, then the asexual mean fitness 
is given by. 



K(t 



Kyy 



-XP 



1/2 



(20) 



where x = P^^^ or {o^f^Y^'^, depending on whether the 
asexual mean fitness is given by Ka.i or Ka.2- 

In this case, we keep the highest-order terms involving 
'-^vv/Kyy in the numerator and denominator of Eq. (15), 
which gives. 



V jcrit 
IP 



(21) 



This equation is interesting, for it says that when ujyy/Kyy 
is large, then sexual replication will outcompete asexual 
replication as long as Hvv/{ip) < 2(l/x^ ~ 1)^- 

Figure 6 shows a plot of {nvy/{'^p))crit as a function of 
^vv/ i^vv for given values of a, /3, and p. 



IV. SEXUAL REPLICATION WITH DISTINCT 
SPERM AND EGG GAMETES 

A. The need for distinct gametes 

As organism size increases, the maturation time in- 
creases as well. If we let N denote the number of cells in 



the adult organism, then it has been suggested that the 
maturation time is proportional to either N^/'^ [13, 13 or 
N^/^ [H, [HOl. To account for both possibilities, we take 
the maturation time to scale as A^", where a = 1/4, 1/3. 
This then implies that Kyy oc iV~". 

If each cell takes up a volume then the total volume 
of an adult organism is Nv. Therefore, if there are n 
adults in the population, then the total volume is on the 
order of nNv, so that the population density is on the 
order of l/{Nv) oc 1/7V. So, p oc 1/iV, which gives that 
Kyy/i'yp) oc N^^". This goes to oo as oo. 

Now, for the sexual strategy to outperform the asexual 
strategy as organism size increases, it is necessary that, as 

A^ ^ oo, {Kvv/i'JP))cnt > Kyy/{jp). Since Kyy/{-/p) 

oo as A^ — > oo, it follows that {kvv / {'yp))crit ^ oo as 
A^ — > oo. Since this can only occur in the small ujyyj Kyy 
region, it follows that we want, 



^™ ^ o (l-x)^ IP _^ ^ o (l -X)^ 

< 2 ^ UJyy < 2 7/9 

Hvv X P l^vv X P 



(22) 



so that, since 7p oc 1/A^, it follows that ojyy should de- 
crease at least as quickly as 1/A^. 

Therefore, for sexual replication to outcompete asex- 
ual replication as organism size increases, the organismal 
sporulation rate should decrease as at least the reciprocal 
of organism size. This is clearly too small compared with 
actual organismal sporulation rates, so that the identical- 
gamete sexual replication model presented in this section 
is insufficient to account for the observation that com- 
plex multicellular life prefers a sexual replication strategy 
(and relies on it exclusively in certain cases). 

Of course, this highly simplified analysis ignores the 
connection between organism size and organism complex- 
ity, and hence the dependence of genome size, and of p 
and x^ on A^. Nevertheless, it is interesting that one does 
not readily obtain from a simplified model that sexual 
replication should be the preferred replication strategy 
as organism size and complexity increases. As will be 
subsequently shown, the assumption of distinct gametes 
provides a simple solution to this problem. 



B. The basic model 

We consider a sexually replicating population of or- 
ganisms that produce two distinct types of gametes. We 
assume that one of the gametes is much smaller than the 
other. The larger gamete, which we term the egg, con- 
tains, in addition to its haploid complement of chromo- 
somes, all of the necessary material to sustain the growing 
embryo during its development from a single cell into an 
immature organism. The smaller gamete, which we term 
the sperm, contains little more than the haploid comple- 
ment of chromosomes. However, the smaller size of the 
sperm means that they can be produced at a far greater 
rate than the eggs. 

We will also assume that a fraction A of the popula- 
tion produces sperm, and a fraction 1 — A produces eggs. 
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at any given time. This assumption does not necessar- 
ily assume a male-female split, since hermaphrodites can 
divide their investment in the production of sperm and 
eggs such that the total output corresponds to an effec- 
tive fraction A of the population producing sperm, and 
the remaining fraction 1 — A of the population producing 
eggs. 

To adapt our notation to this distinct gamete model, 
wc let Us^v, "e,u denote the number of viable sperm and 
egg haploids, respectively. We also let (jJs,vvj ^e,vv denote 
the production rate of diploid spores that then proceed 
to develop into sperm and egg gametes, respectively. All 
other terms are unchanged. 

We also assume, for simplicity, that the sperm cells 
do not have any flagella or any other mechanism to allow 
them to actively transport themselves through the water. 
Rather, we assume that the sperm are dispersed through- 
out the aqueous medium by convection and Brownian 
motion. We will discuss this assumption, as well as a 
number of others, toward the end of the paper. 

Finally, in contrast to the identical-gamete sexual 
replication model presented in the previous section, here 
we will generalize our model somewhat and allow for hap- 
loid death. To this end, wc let k^.^, K^.e denote the first- 
order decay constants of the sperm and egg, respectively. 

With these definitions in hand, the ordinary differen- 
tial equations governing the evolutionary dynamics of the 
population is given by. 



dt 



VV '"Sl.fVV 



(23) 



If we make the same assumptions about organism den- 
sity and work with the same mean fitness and population 
ratios as the previous sexual model, then we have, 



dXai', 



dt 
dt 
dt 



1Xu)s,vv'P ^ pXs,v^e,v {^^d,s ~\~ '^(^))'^s,?; 

2(1 - \)uJe,vvP - lP3Cs,v^e,v " {i^d.e + ^{t))Xe,v 

(24) 



In Appendix B, wc will discuss how this model is af- 
fected if we assume that the egg, sperm, and immature 
organisms themselves take up a certain volume. 

C. Steady-state mean fitness 



Using the fact that Xsi,vv = K{t)/Kvv, we have, at 
steady-state, that. 



IP _ , = Oc) ^ K.(t = .3C) 

Xs,vXe,v — \-L ~r ) 

Kyy f^VV ^VV 



2 \ ^^^^^ — fl 4- ^l(£^22) \ M*^22) 



Xs.V — 



2(1 - A)i^p - (1 -h B(*-°°))B(*='>°) 



The bottom two equations may be plugged into the 
first and re-arranged to give, 



(2A^p-(l 



K{t — Oo) \ K{t — Oo) 



)(2(1-A)^p-(1 + 



K(t — Qo) \ K{t — 00) 



W K.{t = 00) . Kd,, 



K(t= 



^(1 + 



K(t—Oo) 



Kyy 

7P 



(26) 



This equation may be solved to give the steady-state 

mean fitness of the population. In particular, when there 
is no cost for sex, i.e. when Kw/ {jp) 0, then we obtain 
that the steady-state mean fitness must satisfy one of the 



following equations: 

K(f = 00) 2 K{t = OO) ^^ UJs,vv „ 

( ) +{ )-2A p = 

Kyy Kyy Kyy 

Kyy Kyy Kyy 

(27) 

Because in general 

^e,vv "^"^ ^s,vvj we expect that 
Aaig << (1 — yjLUg yy. This implies that the second 
equation yields the smaller steady-state mean fitness, and 
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so this is the steady-state mean fitness of the population 
(essentially, because the eggs are released at a lower rate 
than the sperm, the eggs are the rate-limiting factor con- 
trolling fitness). 
We then obtain, 



fi{t - oo) = ^^^(-1 + ./l + 8(1 - A)^p) (28) 



V. COMPARISON OF THE VARIOUS 
REPLICATION MECHANISMS 

A. Parameter behavior in the limit N oo 

To show how a distinct-gamete sexual replication strat- 
egy allows for sexual replication to be the preferred repli- 
cation strategy in large, complex organisms, we need to 
first explore how the various parameters in our models 
scale with organism size. The parameters we will consider 
are (1) Asexual spore production; (2) Identical-gamete 
spore production; (3) Egg production; (4) Sperm pro- 
duction; (5) The coupling parameter 7. 

We consider each of these parameters in turn: 



1. Asexual spore production 

Complex, multicellular organisms contain genomes 
that encode for a multicellular survival strategy. This 
survival strategy is one that should be more or less opti- 
mized for exploiting a certain environmental niche, oth- 
erwise, natural selection would eliminate it. As a result, 
until the organism reaches its mature, adult size, it is not 
at maximal fitness. Furthermore, because the multicel- 
lular organism is a highly interconnected network of dif- 
ferentiated cells, the organism must contain a minimum 
number of cells before it can even implement its multicel- 
lular survival strategy. Until the organism reaches this 
stage of development, it is not able to exploit the niche 
for which it was designed, nor can it optimally function 
in other niches, since its genome was not designed for 
them. 

To illustrate, a single human fertilized egg, although it 
is an undifferentiated cell, is far less suited to surviving in 
the wild than say, a Paramecium, which was specifically 
designed for this purpose. The reason for this is that 
the fertilized egg contains numerous genes involved in 
implementing a multicellular survival strategy, genes that 
are completely unnecessary for a free-living cell. These 
genes code for various biochemical pathways that cost 
time and energy, and do not contribute directly to the 
production of new cells. As this embryo develops in the 
womb, it becomes even more vulnerable, for it is now in 
a transition from a single cell to a multi-celled organism. 
During this intermediate stage, the developing human 
is neither able to function as a free-living cell, nor as a 
multi-celled human organism. 



Because an organism is vulnerable until it has reached 
a minimal stage of development, complex organisms must 
provide a protected environment for a single-celled spore 
that allows it to develop into an immature organism that 
can then exploit the niche for which it was designed. This 
protected environment must include sufficient quantities 
of pre-processed nutrients to allow the organism to de- 
velop to a certain size and complexity (so that the or- 
ganism does not have to procure resources from niches 
for which it is poorly designed), as well as provide a de- 
fense against harsh environmental conditions, predators, 
and pathogens. In Appendix C we provide a simplified 
analysis of these assumptions. 

For many organisms, the way to provide such a protec- 
tive, nutrient-filled environment is to produce relatively 
large eggs that contain an outer protective shell, as well 
as a "yolk" that provides sufficient nutrients for the ini- 
tial, single-celled spore to develop into an immature, yet 
fully differentiated organism that can exploit the envi- 
ronmental niche for which its genome was designed. 

In models for the growth of an organism to full size, 
which have shown a power-law dependence of growth 
time with the number of cells in the organism, the un- 
derlying assumption is that the metabolic rate of the or- 
ganism is diffusion-limited. The 1/3 power-law comes 
from assuming a spherical model, while the 1/4 power- 
law comes from assuming a fractal geometry to the organ- 
ism, consisting of many branching blood vessels spread- 
ing throughout the body [13, 113 • 

In any event, assuming that an asexually replicating 
adult organism is fully devoted to producing eggs, then 
if / denotes the ratio of the number of cells in the im- 
mature organism to the number of cells in the adult, we 
obtain that egg production may be modeled as a pro- 
cess whereby an A^-celled adult grows to (1 -I- f)N cells, 
and then the material from the fN cells are discharged 
to produce an egg (in reality, these fN cells consist of 
the "yolk" material necessary to produce the fN - celled 
organism later on). 

Assuming that the rate of growth of the adult is pro- 
portional to the metabolic rate, with a constant of pro- 
portionality 77, that the metabolic rate scales as N-^^", 
and that / is small, we obtain, 

regg = (29) 

where r^gg denotes the characteristic time it takes to pro- 
duce a single egg. 

If we let uJae,vv denote the production rate of diploid 
eggs in the asexual population, then we have that 
u;ae,vv oc l/r^gg, so that, 

^ae,vv CX TV"" (30) 

It should be noted that there is a class of organisms, 
which are termed viviparous, that do not produce rela- 
tively large eggs that are released into the external envi- 
ronment, but rather give birth to live offspring. In these 
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organisms, the eggs are much smaller, containing httle 
to no yolk. Instead, the organisms develop inside the fe- 
male for a certain period, and then emerge as immature 
organisms that continue to develop on their own. 

Although this paper considers broadcast fertilization, 
we argue that the existence of viviparous organisms does 
not contradict our basic argument for the existence of 
eggs. In viviparous organisms, the female can afford to 
produce eggs with little to no yolk, because she uses her 
body as the protective environment in which her offspring 
can develop. Furthermore, she directly provides her off- 
spring with the nutrients necessary to grow and develop. 
Therefore, the total required investment of resources in 
viviparous organisms is similar to that of oviparous (egg- 
laying) organisms (indeed, it is likely that the resource 
cost in viviparous organisms is higher, since the female 
must carry the offspring inside of her, which limits her 
mobility). As a result, both viviparous and oviparous 
organisms have similar constraints on the rates at which 
females can produce offspring, leading to similar scaling 
of offspring production rates with size. 

We should note that, if we consider parental care, then 
this may be regarded as an additional input of resource 
that further reduces the effective reproduction rate of 
females. 



2. Identical- gamete spore production 

With the identical-gamete sexual replication model, 
the organisms produce haploid eggs which must then re- 
combine with one another to form a diploid, "fertilized" 
egg that then develops into an immature organism. 

We note first that if a diploid egg contains enough ma- 
terial to produce fN cells in the immature adult, then 
each haploid egg contains enough material to produce 
fN/2 cells. 

There are two ways to produce the haploid eggs: (1) 
The full-sized diploid egg divides in two, in which case 
the value of ujvy for this model, which we denote ujie,vv, 
is simply ujae,vv (where "i" stands for "identical"); (2) 
An initial diploid cell divides into two haploids, and each 
haploid is sequentially filled with enough material to form 
the two haploid eggs. When / is small, then the previous 
analysis also yields that LUi^^vv — ^ae,vv 

Therefore, in any event we may assume that oJie^w = 



3. Egg production 

With the distinct-gamete sexual replication model, the 
egg-producing organisms produce haploid eggs that are 
the same size as the diploid eggs in the asexual organisms, 
since the sperm contain essentially only genetic material, 
and no additional material to supply the growth and de- 
velopment of the embryo. 



Following a similar analysis to the one carried out 
above, we may then show that, if / is taken to be small, 
then LiJe^vv — ^ae,t)i;/2. The factor of two comes from the 
fact that the haploid egg in this case must contain twice 
the material of the haploid eggs in the identical-gamete 
case, which reduces the egg production rate by a factor 
of two. 



4- Sperm production rate 

Since sperm do not contain any "yolk" material, then, 
assuming that an adult sperm-producing organism is 
fully devoted to producing sperm, we may take oJs.w 
to be proportional to the metabolic rate of the organ- 
ism. Therefore, we may assume that liJs,vv scales as 

For generality, however, we take uJs,vv oc 
where (3 G (0, 1). 

5. The coupling constant 7 

The second-order rate constant that determines the 
haploid fusion rate depends on the size of the haploids. 
On the one hand, larger haploids have a larger sur- 
face area available for coupling. On the other hand, 
larger haploids will have a slower rate of random motion 
through the aqueous environment, which reduces the fu- 
sion rate. 

Following the approach of Dusenbery [2^, we model 
haploid fusion as a binary collision process between ran- 
domly moving particles in a medium. Since in this paper 
we do not assume that the sperm have any kind of active 
transport mechanism (e.g. flagella), the random motion 
is due to Brownian motion. 

Based on the theory of chemical reaction kinetics, it is 
possible to show that, 

7 = 7r(ri-f r2)2— ^— ^ (31) 

where ri , vi denote the radius and average speed of the 
small particles, and r2, V2 denote the radius and average 
speed of the large particles. 

In the identical-gamete model, ri = r2 = re, and vi = 
V2 — Ve, where r^. and are the radius and average 
velocity of the eggs. This gives. 



where "ig" is meant to signify identical gamete. 

Now, the radius of the egg, re, is related to its volume 
via Ve = (4/3)7rrg, and the egg volume is proportional to 
fN/2. Therefore, oc N^/^ ^ rl (x N^/^. 

Regarding Ve , we note that diffusive motion of the eggs 
through the aqueous medium occurs because of random 
collisions of water molecules against the egg surface. This 
leads to a random stochastic force F{t), with an average 
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magnitude F. If we assume that velocities are sufficiently 
small so that we are in the laminar-flow regime, then for 
a spherical egg the drag depends on velocity via D = 
QiTfireVe, where ^ is the fluid viscosity [26| . Equating F 
and D gives a velocity, 



(33) 



and so Ve is proportional to l/^e, and hence Ve (x N 
Putting everything together, we obtain that. 



(34) 



In the distinct-gamete model, we have ri = Vg, r2 — re, 
where is the radius of the sperm, and is the radius 
of the eggs. We also have vi = Vg, and V2 = Ve- Since 
the eggs are much larger than the sperm, we have « 0, 
v„ ~ 0, so that. 



(35) 



where "dg" is meant to signify distinct gamete. 

Since the sperm remain more or less of constant size 
as organism size increases, the only quantity that scales 
with N is re, giving, 



7dg OC 



(36) 



B. 



How parameter behavior affects the mean 
fitness as A'^ ^ oo 



1. Identical gametes 

In the identical gamete model, we have that Kvv/{^p) 
scales as A^^/s-a ^ cx) as iV ^ oo, since a — 1/4,1/3. 
This implies that LUie,vv must decrease at least as quickly 
as 7p OC iV^^/'^. Since tOie^w — ^ae,vv oc N^", where 
a = 1/4, 1/3, it appears that the actual decrease in egg 
production rates as a function of organism size is not 
sufHciently rapid for the sexual replication strategy to 
outcompete the asexual one. 



2. Distinct gametes 

In the distinct gamete model, we have that Km/iip) 
scales as iV^/^~". Therefore, Kw/i'jp) either remains 
essentially constant as — > oo, or scales as iV^/^^. How- 
ever, because we have distinct sporulation rates (jJs,vv and 
uje,vv, we cannot directly apply the previous analysis to 
this model. 

To determine how the mean fltness scales as a function 



of organism size, we flrst re-write Eq. (26) as, 

2(1 _ X)'^p - (1 + ~<t-°°)fit-°°) 

1 1 K(t — Oo),^ K(t — oo) ^ 

^ -{! + — -) 



X {fi{t ^ oo) + Kd,s){K{t = oo) 4- Kd,e) 
^ 



(37) 



The normalized mean fitness K(t = oo)/ Kw is a de- 
creasing function of Kyy/{'jp), and so has a maximal value 
given by Eq. (28). We therefore have that, 

ii(t — oo) , Rit — oo) , , , „„ 
^(1 + '-) < 2(1 - X)-^p 



K.{t — oo) < 2(1 — X)LJe.vvP 



< 



Aw, 



2X^^p -(1+ - XuJs,vv - (1 - A)We,™ 

(38) 

These inequalities may be plugged into Eq. (37) to give, 


< 2(1 - A)^p - (1 + -<^ = ^)f^^ = ^) 

< — (i _ 1) 

X (2(1 - X)LL!e,vvP + '«<i,s)(2(l - X)LUe,vvP + Kd,e) 

^ (39) 



l„(i 



Now, as iV — > oo, uJe,vv and scale as A'' Ws,vv 
scales as , 7 scales as N^^^, and p scales as . This 
implies that 1/(1 — (1/A — l){uje.vv/i^s,vv)) — > 1 as ^ 
00. Furthermore, if we assume that Kd,s,fie,s > 0, then 

(2(1 - \)uJe,vvP + K(i,s)(2(l - \)uJe,vvP + Ke^s) ^ Kd,sKe,s 

'AS N 00. 

Therefore, as A^ — > 00, the right-most term in Eq. (39) 
scales as N^^^~^, which goes to as long as /3 > 1/3. 
In this case, the cost for sex disappears as organism size 
increases, so that the steady-state mean fitness is given 
by the mean fitness when there is no cost for sex. 

The assumption of /? > 1/3 certainly holds if /3 = 1 — a, 
for then we obtain that (3 = 2/3,3/4 > 1/3. This lower 
bound seems to be far below the actual rate at which the 
fraction of sperm-producing cells in an organism scales 
with organism size. To illustrate with a simple order of 
magnitude calculation, the number of sperm in the aver- 
age ejaculate of an adult human male is approximately 
2.5 X 10® while the number of sperm in the aver- 
age ejaculate of an elephant is approximately 2 x 10^^ 
[52| . If an average adult human male has a mass of 
100 kg, while an average adult elephant male has a mass 
of 5,000 kg [s^l, then the average ejaculate of an adult 
elephant should contain 1.25 x 10^° sperm. 
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Therefore, the adult elephant seems to produce an or- 
der of magnitude more sperm than what would be ex- 
pected assuming only a linear scaling law, suggesting that 
/3 > 1, at least for this example. While we have not taken 
into account issues such as frequency of ejaculation, this 
result nevertheless suggests that the requirement that 
/3 > 1/3 is easily satisfied for complex organisms. 



C. Comparison of the three replication mechanisms 

As TV cx), the identical-gamete sexual replication 
strategy becomes too costly for it to remain competi- 
tive with either the asexual or the distinct-gamete sexual 
replication strategies. However, under the very loose re- 
quirement that P > 1/3, the cost for sex in the distinct- 
gamete strategy decreases to zero as oo, so that the 
steady-state mean fitness of the distinct-gamete strategy 
for large N is given by. 



^^4(-l + ./l + 8(l-A)^^p) (40) 

The elimination of the cost for sex is due to the sperm 
gametes. Because they are small and can be produced in 
large quantities compared to the eggs, for large organisms 
the eggs essentially become envelo ped in a "sperm cloud" 
that ensures rapid fertilization [1. l22j|. 

Using the fact that uje^y^ = L0ae,vv/'2, we have that the 
asexual mean fitness is given by. 



Ka{t = Oo) 1 

= max| -(— 1 - 

hhyy 2 



-'ae,vv 2 



p2), 



4a/3- 



-P)} 
(41) 



while the sexual mean fitness is given by. 



~ = ^(-1 + Jl + 4(1 - A)^^p) (42) 

Note then that if p < 1 — A, the sexual mean fitness 
outcompetes the asexual mean fitness when the asexual 
mean fitness is given by the first fitness function. Also, 
if /3 < 1 — A, then the sexual mean fitness outcompetes 
the asexual mean fitness when the asexual mean fitness 
is given by the second fitness function. 

Therefore, given values for p and /3 that are less than 1, 
we can choose a A £ (0, 1 — max{p, /?}), and obtain that 
sexual replication will outcompete asexual replication for 
sufficiently large organisms. Of course, the smaller the 
required value of A, the fewer the number of organisms 
that are involved in producing sperm, and so the larger A^ 
has to be before there is a sufficient amount of sperm to 
reduce the time cost for sex to an extent that it becomes 
advantageous over the asexual strategy. 

For the case where there is an equal number of sperm 
and egg-producing organisms, sexual replication will only 



outcompete asexual replication for large organisms if 
p, /3 < 1/2. This suggests that a sexual rephcation 
strategy employing distinct sperm and egg gametes will 
outcompete an asexual strategy for large organisms if 
the replication fidelity is sufficiently low, and if the fit- 
ness penalty for having a partially defective genome is 
sufficiently large. Given that genome length increases 
with organismal complexity, and given that the genome 
encodes for an interconnected network of biochemicals, 
cells, and organs that are crucial for organismal survival, 
it is likely that the conditions given above hold as A^ — *■ oo 
(though admittedly, we do not have actual data to verify 
this statement). 



VI. DISCUSSION 

A. The evolution of distinct gametes, distinct 
sexes, and the sex ratio 

The results of this paper suggest that a distinct-gamete 
sexual replication strategy leads to the selection for sex- 
ual replication over asexual replication in larger organ- 
isms. Furthermore, if the replication fidelity is suffi- 
ciently low, and if the fitness penalty for having a par- 
tially defective genome is sufficiently high, then the fit- 
ness benefit for sexual replication can even be shown to 
overcome the two-fold cost for sex. 

In the context of sexual replication via broadcast 
fertilization, distinct sperm gametes may have evolved 
as "parasitic" gametes defecting from an egg-producing 
strategy. The idea is that, in a population consisting en- 
tirely of egg-producing organisms, an organism that pro- 
duces smaller gametes in greater numbers will have a se- 
lective advantage, since it can fertilize many eggs for the 
same cost of producing one egg. This results in an evolu- 
tionary pressure for evolving small gametes that can be 
produced in great numbers, followed by a co-evolutionary 
pressure for evolving large gametes that contain the nec- 
essary "yolk" material for allowing the fertilized egg to 
develop into an immature organism. 

Broadcast fertilization is the earliest sexual replica- 
tion strategy in complex organisms. This potentially ex- 
plains why organisms that do not replicate via broad- 
cast fertilization nevertheless employ distinct sperm and 
egg gametes: The distinct-gamete strategy proved advan- 
tageous for this type of fertilization mechanism, which 
essentially "locked" this strategy into subsequent evolu- 
tionary lines. While it is in principle possible for an or- 
ganismal line to have evolved directly from asexual repli- 
cators to sexual replicators via mating, it may be that 
the body plan and behaviors necessary for mating re- 
quire an organism more complex than one that sexually 
replicates via broadcast fertilization. As a result, it may 
be far more likely that an organismal line that evolves 
from asexual replicators to sexual replicators via mating, 
does so via an intermediate that replicates sexually via 
broadcast fertilization. 
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Alternatively, it should be noted that the same evolu- 
tionary pressures driving the formation of distinct sperm 
and egg gametes for broadcast fertilizers also apply to 
mating organisms. Therefore, there may be a selective 
pressure for anisogamy even for a sexual replication strat- 
egy that involves mating. The key difference between the 
selective pressure for anisogamy in the case of broadcast 
fertilization, versus the selective pressure for anisogamy 
in the case of mating, is that in the former case anisogamy 
is necessary for maintaining a selective advantage for the 
sexual strategy itself. In the case of mating, where the 
time cost associated with haploid fusion is negligible, 
anisogamy outcompetes other sexual replication strate- 
gies, but is not necessary for maintaining a selective ad- 
vantage over asexual replication. 

With distinct sperm and egg gametes, it is possible to 
have distinct sexes, where one sex, the male, produces 
only sperm, while another sex, the female, produces only 
eggs. The advantage of such a strategy is that since each 
sex can focus on producing only one type of gamete, the 
overall efficiency of gamete production is increased. The 
disadvantage of such a strategy is that replication is de- 
pendent on the mating of a male and a female, and so can 
have a large fitness penalty if the contact frequency be- 
tween males and females is too low. In such cases, it may 
be advantageous for the organisms to be hermaphrodites, 
and to be capable of engaging in self-fertilization when 
necessary. Indeed, the disadvantage in having distinct 
sexes has led to many organismal lines preserving a min- 
imal ability for asexual replication, when necessary (it 
has recently been discovered that female sharks can re- 
produce asexually when not in contact with any males for 
a sufficiently long period of time). In other organismal 
lines, the organisms can change sexes so as to maintain an 
appropriate ratio of males to females to make the sexual 
strategy the preferred one. 

It is now established that mammals are the only class 
of organisms known to be obligately sexual, with distinct 
male and female sexes jl3<] . Presumably, mammals are 
sufficiently complex that the costs associated with sex- 
ual replication are outweighed by its benefits, to an ex- 
tent that maintaining both sexual and asexual replication 
pathways, as well as an ability to switch between them, 
simply incurs a fitness penalty due to the costs involved. 

A major issue in evolutionary biology is to understand 
why the sex ratio in many populations is close to 1 : 1. 
The generally accepted theory holds that the sex ratio is 
approximately 1 : 1, because it is an Evolutionarily Sta- 
ble Strategy [5^. The basic argument is that, although 
there are many organisms in which most of the males 
of the population never mate, it is nevertheless advan- 
tageous for half of an organism's offspring to be male, 
since the few male offspring that do mate will mate with 
many females and father numerous children. Further- 
more, the production of excess males allows for compe- 
tition between individual males to claim a group of fe- 
males. This competition, if not too costly, provides a 
natural mechanism to select for beneficial genes that are 



then passed on to offspring. 

While this paper does not focus on the issue of the 
sex ratio, we nevertheless wish to bring up two possible 
alternative explanations for the 1 : 1 sex ratio. In doing 
so, we should first note that one has to be careful in 
stating that a 1 : 1 sex ratio is an ESS, since in certain 
kinds of organisms, such as ants and bees, the sex ratio 
is far from 1 : 1 (a colony consists of a single queen, 
numerous sterile female workers, and few male drones 
that mate with the queen and fertilize her eggs). 

Nevertheless, in certain contexts, a 1 : 1 sex ratio may 
prove advantageous. For example, in an environment 
where the population density is low, a 1 : 1 sex ratio 
maximizes the contact rate between males and females, 
and therefore minimizes the fitness penalty associated 
with a distinct-sex strategy. To see this, note that, in a 
population of N organisms where the fraction of males is 
X, the number of males is xN and the number of females 
is (1 — x)N, so that the contact rate between males and 
females is proportional to x{l — x)N'^ , a quantity that is 
maximized when x = 1/2. Given that larger organisms 
typically have smaller population sizes, it may be that 
species that maintain an approximately 1 : 1 sex ratio 
do so because they evolved along organismal lines with 
relatively small populations. 

It is also possible that the 1 : 1 sex ratio is a conse- 
quence of the fact that certain organisms have specific 
sex-determining chromosomes, e.g. "X" and "Y" . Since 
each organism is produced by one sperm and one egg, on 
average one would expect half of all offspring produced 
to be male, and half female. 

The use of specific sex-determining chromosomes may 
not necessarily lead to an optimal sex ratio. However, 
this method for sex determination may be a relatively 
simple one in certain environments, so that it might be 
more robust and less costly than other methods. As a 
result, while other sex-determining strategies may lead 
to more optimized sex ratios, the additional costs they 
entail exceed the fitness benefit that they provide, so that 
such strategies are at a net disadvantage. 



B. Group selection versus individual selection 

In this paper we have relied on the use of the mean 
fitness of a population to determine the replication strat- 
egy that is advantageous in a given parameter regime. 
This approach is known as the group selection approach, 
because it assumes that genes that are beneficial for the 
population as a whole are the ones that will be selected. 

Strictly speaking, because an individual organism is 
the reproducing agent, it is not necessarily true that a 
strategy that benefits the group is the one that benefits 
the individual. Indeed, one of the central points of game 
theory is that individuals acting in their own self-interest 
can often engage in behaviors that are detrimental to the 
group as a whole. The Prisoner's Dilemma is an excellent 
example of this. Biologically relevant examples include 
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the emergence of cancer in multicellular organisms, vi- 
ral evolution, and possibly even the formation of neural 
pathways that lead to addictive behaviors [i^ . 

In general, if the individual organisms in a group are 
competing for limited resources, so that one individual 
increases its fitness at the expense of the other organ- 
isms, then the group selection approach will not correctly 
predict what strategies will be selected. However, if in- 
dividual organisms can maximize their fitness without 
adversely affecting the fitness of other organisms, then 
the group selection approach will correctly predict what 
strategies will be evolutionarily selected. 

In the case of the sexual replication models being con- 
sidered in this paper, the group selection approach is 
expected to be valid, because the production of sperm, 
although it may arise as a form of defection from an in- 
vestment in egg-production, nevertheless serves to reduce 
the time cost for sex. As a result, an organism that pro- 
duces sperm gametes is not necessarily increasing its own 
fitness at the expense of other organisms. On the con- 
trary, a female that invests in the production of haploid 
eggs that are fertilized by the sperm pays a small time 
cost for sex. By contrast, a female that produces hap- 
loid eggs that fuse with other haploid eggs pays a large 
time cost for sex. As a result, although the former fe- 
male produces eggs half as quickly as the latter female, 
the former female's eggs are fertilized much more quickly, 
and therefore her overall reproduction rate is higher. 

Of course, this argument assumes that sperm competi- 
tion is negligible, or that the fitness benefit of sperm com- 
petition outweighs the cost. Once this no longer holds, 
then sperm competition will lead to a co-evolutionary 
dynamics that will result in a reduction of the overall 
population's mean fitness. However, as long as this re- 
duction is not so large so as to make asexual replication 
the preferred replication strategy, the sexual replication 
strategy will still dominate. 



VII. CONCLUSIONS AND FUTURE 
RESEARCH 

This paper analyzed the evolutionary dynamics asso- 
ciated with three distinct replication strategies for mul- 
ticellular, sporulating organisms: (1) A purely asexual 
strategy; (2) A sexual strategy employing identical ga- 
metes; (3) A sexual strategy employing distinct sperm 
and egg gametes. Under the assumption that the sex- 
ual populations replicate via broadcast fertilization, we 
found that the distinct-gamete strategy is necessary for 
maintaining a preference for sexual replication as organ- 
ism size increases. 

As was mentioned earlier in this paper, previous stud- 
ies exploring the selective advantage for a distinct-gamete 
strategy have focused on the preference for a distinct- 
gamete sexual replication strategy over an identical- 
gamete sexual replication strategy. This paper, by con- 
trast, argues that a distinct-gamete strategy is necessary 



for maintaining the selective advantage for the sexual 
replication strategy itself. We believe that this is an 
important result, for a theory explaining the preference 
for one sexual replication strategy over another does not 
explain why sexual replication should exist in the first 
place. 

The analysis of this paper relies on the assumption that 
complex organisms must replicate by producing compar- 
atively large eggs. As a result, the necessity for produc- 
ing eggs cannot be dependent on the selective advantage 
for the distinct-gamete sexual replication strategy itself. 
Otherwise, if large eggs were not necessary for organis- 
mal viability, then it would make much more sense for 
organisms to replicate by releasing microscopic single- 
celled spores. However, according to the analysis in [2l|, 
in this case the asexual pathway would be the advanta- 
geous one, and so there would be no selective pressure 
for evolving any kind of sexual replication strategy. 

The models considered in this paper are highly sim- 
plified, and so future research should involve developing 
more realistic replication models. Some of the features 
that should be considered are: (1) More realistic fitness 
landscapes, derived from organismal genomes consisting 
of multiple genes and more than two chromosomes. Fur- 
thermore, the assumption that a single mutation renders 
a single gene or genome region defective needs to be more 
closely examined. (2) Gamete-release cycles, whereby the 
organisms do not release their gametes continuously, but 
rather store up their gametes for a certain period and 
then release them rapidly and in large quantities. Pre- 
sumably, such a strategy can lead to high gamete den- 
sities, resulting in high gamete contact rates and rapid 
fertilization. This in turn could lead to a more rapid de- 
crease in the cost for sex with increasing organism size 
than our model predicts. (3) Sexual replication models 
involving explicit mating between organisms, either via 
external fertilization or internal fertilization. (4) Mobile 
sperm gametes, and eggs that produce pheromones to 
attract sperm. (5) Mortality. 

Finally, part of the difficulty in understanding the evo- 
lution and maintenance of sex is that there is a wide 
variety of sexual and mixed sexual-asexual strategies em- 
ployed by different species, making it difficult to develop 
a single model for analyzing the selective advantage for 
sex that is applicable to all organisms. By studying the 
time and energy costs, as well as the fitness benefits, 
of the replication strategies employed by different organ- 
isms, it may be possible to obtain a better understanding 
both of the conditions under which sexual replication is 
an advantageous strategy, as well as the reasons why or- 
ganisms adopt the replication strategies that they do. 
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APPENDIX A: THE DEPENDENCE OF THE 
MEAN FITNESS K{t = oo) ON p IN THE 
ASEXUAL REPLICATION MODEL 

We wish to prove that there exists a Pcrit G P] 
such that Ka,i > Ka^2 for p e {pcrit, 1] and Ko,i < '«a,2 for 
p G (0, Pcrit)- We will do this in three steps: 

1. We will prove that > i^a,2 for P S (/3, 1]. 

2. We will prove that Ka,i < Ha,2 for p € (0, af3). 

3. Wc will prove that there exists a unique pcrit S 
[aj3,j3] such that Ka,i > Ka,2 for p G (Pcriti/?] and 

Ka,l < Ko,2 for P e [a/3, Pcrit)- 

Note that all three statements together complete the 
proof. 

To prove the first step, consider the function /(.x, y) = 
—X + v'a;^ + xy. Note that / is an increasing function of 
y for x,y > 0. Note also that, 



9/ _ x + I - y^sF+xy 
dx A/a;2 + xy 



(Al) 



Now, for x,y > 0, {x + yjl) = x + xy + y /4 > 
x^ + xy =^ x + y/2 > ^Jx'^ + xy, so that df /dx > for 
x,y > 0, and hence / is also an increasing function of x 
for x,y > 0. 

Therefore, given xi,X2 such that < xi < X2 and 
yi,y2 such that < yi < 1/2, we have f{xi,yi) < 

If we set xi = a, X2 = 1, yi = 4:{LUyy/Kyy)(3p, 2/2 = 
4:{ujyy / Kyy)p'^ , thcu ioT p > (3 wc havc < a;i < a;2 and 
< 2/1 < 2/2, and so. 



— = -(-1 + ./1 + 4— p2) 



> ^(-a + Wa' + 4— a/3p) = ^ 

KaA > Ka,2 



(A2) 



thereby proving the first step. 

To prove the second step, we will prove the equivalent 
statement that Ka,il Kyy — K,a,2/i^vv < for p £ (0, a/3). 

To this end, define, for given p, a, (3, the function g{X), 
via, 



g{X) = -(1 - a) + + Ap2 - v'a2 + Xa^p (A3) 

and note that Ka. i/i^w — >^a.2l*^vv = (1/2).9(A = 
4:{cOyy / Kyy)) . Thcrcfore, if we can show that g{X) < 
for p e (0, a/3), then this will establish that Ka,i/Kyy — 

Ka, 2/^1,1, < for p G (0,a,/3). 

Now, note that g(0) = and that. 



dg _ I ,P\J (x^ + Aa/3p — a/J-^l + Ap^ 
^^A ~ 2^^ v'(l + Ap2)(a2 + Aa/3p) 



We claim that dg/dX < for p e (0, a/3). To see this, 
note that. 



<^ pV + Aa/3p < a/3-v/ 1 + Ap^ 
^ p^o? + Aa/3p3 < a^^^ + Xo?0^p^ 
^ o?{0^ - p^) + Xaf3p'^{af3 - p) > 



(A5) 



This last statement holds for p e (0,a/3), since this im- 
plies that < p < a/3 < /3. 

Therefore, the identity g{0) = and the inequality 
g'{X) < for p G (0, a/3) together imply that 5(A) < 
for p G (0, a/3) and A > 0, and so the second step has 
been proved. 

By continuity, the first and second statements imply 
that Ka,i > Ka,2 for p = /3, and Ka,i < Ka,2 for p = a/3. 
To prove the third and final step, we define A = 

'i:{Wyy / Kyy) , SO that, 

Ko,l = Ka,2 

<^ -1 + + Ap^ = -a + ^/a'^ + Xa(}p 
a/1 + Ap2 - -\/a2 + Aa/3p = 1 - a 



^ Ap2 + Aa/3p + 2a = 2v/(l + Ap2)(a2 + Aa/3p) 
<^ Ap(p - aPf - 4a(l - a)(/3 - p) = (A6) 

Now, when p = a/3, then Ap(p — a/3)^ — 4a(l — a)(/3 — 
p) = — 4a/3(l — a)^ < 0, while when p = /3, then Ap(p — 
a/3)2 - 4a(l - a)(/3 - p) = A/33(l - a)^ > 0. 

Therefore, by the Intermediate Value Theorem, there 
exists a pcrit G [a/3, /3] such that Ap(p — a/3)^ — 4a(l — 

a){(3 - p) = ^ Ka,l = Ka,2- 

To show that Pcrit is unique, we differentiate Ap(p — 
a/3)2 - 4a(l - a) (/3 - p) to obtain A(p - ajSf + 2Xp{p - 
a/3) + 4a(l - a) > for p > a/3. 

Therefore, Ap(p — a/3)^ — 4a(l — a)(/3 — p) is an in- 
creasing function for p > a/3, and hence pcrit must be 
unique. 

Now, if Ka,i < Ka,2 for some p G {Pcrit,P], then 
since Ka,i > Ka,2 for p = /3, it follows that there ex- 
ists some p G (pcrit,/?] such that Ka.i = Ka,2, Contradict- 
ing the uniqueness of Pcrit- Therefore, Ka,i > Ka,2 for 

P G (Pcrtt,f3]. 

Similarly, we can show that Ka.i < Ha,2 for p G 
[aj3,pcrit)j completing the proof of the third statement. 

With all three statements proved, our claim is estab- 
lished. 



APPENDIX B: THE EFFECT ON THE COST 
FOR SEX IF THE VOLUME OF EGGS, SPERM, 
AND IMMATURE ORGANISMS IS 
CONSIDERED 

If the volume of eggs, sperm, and immature organ- 
isms is taken into account, then the system volume will 



17 



increase at a greater rate than predicted by the mod- 
els considered in this paper, which assumed that only 
the adult organisms dictated the total volume of the sys- 
tem. This will lead to an increased cost for sex than 
predicted in this paper. As a result, the asexual strat- 
egy will still outcompete the identical-gamete strategy as 
organism size increases. 

However, because the distinct-gamete strategy out- 
competes the asexual strategy as organism size increases, 
it is possible that the increased cost for sex associated 
with egg, sperm, and immature organism volume will 
change this result. 

If each cell has a volume i^, then the total volume oc- 
cupied by the sexually replicating population is, 

V = {nam,vvN + fiN + ne,vfeN + n,,„)j/ (Bl) 

so that, 



rapidly than N^^", so that Xs,v/N scales less rapidly 
than iv^+"-^ If /3 = 1 - a, then x^^y/N scales less 
rapidly than N'^. 

Therefore, as A'^ — > oo, the steady-state solution when 
there is no cost for sex is such that the factor l+fiXai.vv + 
feXe,v + Xs,v/N docs not scale more rapidly than iV°, 
so that Kyy/{'yp) scales no more rapidly with N as it 
does in the model considered in this paper. Since this 
scaling is such that the steady-state solution as ^ cxd 
is the one when there is no cost for sex, we obtain that 
the steady-state solution when there is no cost for sex is 
the self-consistent one as A'' ^ oo. As a result, even if 
we consider egg, sperm, and immature organism volume, 
the steady-state solution in our distinct-gamete model is 
unchanged, and so the conclusions we have drawn from 
the models considered in this paper are unchanged as 
well. 



V 



Nv 



where fi denotes the average ratio of the number of cells 
in the immature organisms to the number of cells in the 
adult, and /e denotes the ratio of the number of cells 
that the eggs contain enough material to produce to the 
number of cells in the adult. 
It is readily shown that, 

~ 1^am,vv (-1 j_ f j_ f I ^ \ 

— p(l "1" fi^ai,vv ~i~ fe^e,v 'J^'^s^v^ {■^'^) 



SO that. 



(B4) 



As a result, if we re- work the distinct-gamete equa- 
tions, we obtain an identical set of steady-state equa- 
tions, except that the cost for sex becomes (k,,,,/(7p))(1 + 
fiXai,vv + feXe,v + {xs,v/N)). This is stiU Kyv/{"fp), except 
that now p is dependent on the steady-state population 
ratios. 

If the steady-state solution is given by the solution 
when there is no cost for sex, then the following equalities 
and inequalities hold: 

Xai,vv ~ ^ 2(1 A) ^ p 

X^^y — 

2\LJs^yyP 2\U}s,VvP 



K{t = OO) -I- Kd,s H{t = OO) 

iXiVs^yyP 



^yy{-l+^l+4^p) 



(B5) 



APPENDIX C: A SIMPLE MATHEMATICAL 
MODEL TO JUSTIFY THE NECESSITY FOR 
EGGS 

In this section, we develop a simple mathematical 
model that illustrates the necessity for egg production 
in complex organisms. We assume that, until a multi- 
cellular organism has reached a minimal size consisting 
of fN cells, where N is the total number of cells in the 
adult and / € (0, 1) is a size fraction, the organism has 
a death rate characterized by a first-order constant Kd- 
While this assumption is overly simplistic, it reflects the 
fact that an organism employing a differentiated, multi- 
cellular survival strategy must reach a minimal level of 
development in order to properly function. 

If the adult organism releases spores with enough ma- 
terial to form an organism with f'N cells, where < 
./' ^ then, following the arguments developed earlier 
in this paper, we have that the time it takes to produce 
a single spore is given by, 



f 

T = —N° 
' spore — ^ 

V 



(CI) 



so that, if LOyy{f') denotes the spore production rate when 
the spore contains enough material to produce f'N cells, 
then u>yy{f) = {'n/f)N-". 

Now, the initial spore grows from an organism with 
f'N cells to an organism with fN cells. To determine 
the characteristic growth time, denoted Tgrow, we note 
that if n denotes the number of cells in the organism, 
then we have, 



dn 



= rjn 



l-a 



1 /" - /" 



' grow 



-N° 



(C2) 



Since the organisms decay with a first-order rate con- 
stant Kd until they reach a size of fN cells, and since 

it takes the newly released spores a time Tgrow to reach 
scales less rapidly than N°, and Xs,y scales less the size of fN cells, the fraction of newly released spores 



As iV OO, uje_yy and Kyy both scale as A^ so 
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that reach a size of fN cells is given by e~'^"'^<'^°^ , so 
that the effective production rate of immature organisms 
is given by, 

^™,e//(/',/) = ^iV-"e-^^^^ 

(C3) 

Note that ujvv.effif, f)/^vvif) = 
{f/r)eM~{^D/v){ir - /'")/a)iV"], a quantity 
that goes to as ^ oo. 

Therefore, this model suggests that as organism size 
increases, it makes sense to produce relatively large eggs 
to maximize fitness. Of course, the arguments presented 



above could be used to justify producing eggs that con- 
tain enough material to produce the full-sized organism. 
What this model neglects is the cost to the fitness of the 
parent by investing heavily in a given offspring. Taking 
such considerations into account would lead to an opti- 
mal egg size that is at some intermediate size. 

Despite the simplicity of our model, we nevertheless 
believe that it illustrates the basic reason for a replica- 
tion strategy based on the production of eggs in complex 
organisms (and for parental care in the most complex 
organisms) . 
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